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I. INTRODUCTION 

Entanglement, a key concept in quantum-information sci- 
ence 0]], has also been playing a pivotal role in the study of 
quantum (or, generally, nonclassical) correlations and related 
consequences in many -body systems Q], especially quantum 
phase transitions (QPTs) induced by varying parameters of 
system Hamiltonian at zero temperature fit] (see also Refs. 
Jj] for geometric approaches to QPTs.). Accordingly, many 
efforts have been made to find the possible relation between 
QPT in a system and nonanalyticities of, e.g., pairwise entan- 
glement of its ground state (GS) JH]. For example, under some 
fairly general conditions, a classification of QPTs in terms of 
pairwise entanglement has been provided in Ref. |6J]. Besides, 
through density functional theoretical arguments, it has been 
shown that by an appropriate choice of entanglement measure 
(either bi- or multipartite), one can study physical properties 
of a system near its critical points 13]. It is thus evident that 
having a sufficiently strong entanglement measure which is 
able to detect entanglement regions of the GS (in terms of the 
Hamiltonian parameters) in a more subtle way is of paramount 
importance. The difficulty of the issue, i.e., quantifying quan- 
tum correlations, is mostly due to the existence of various 
types of entanglement in quantum many-body systems, such 
as free entanglement (FE), bound entanglement (BE) (H], and 
pairwise or multipartite entanglement 121, featuring also new 
phases for quantum matter (e.g., topological order flOlO . Here, 
however, we shall concentrate only on BE. 

Simply speaking, bound entangled states, in contrast to free 
entangled states, are those states which are not directly useful 
for quantum-information tasks such as teleportation Hlltl . Be- 
cause of the positive partial transpose (PPT) property of bound 
entangled states (in the bipartite case) HI], their entanglement 
cannot be distilled by any local operation and classical com- 
munication; for a general review, see, e.g., Ref. jl2]. The 
possibility of the existence of non-PPT BE, however, is still 
an open question lfl3tl . Nevertheless, bound entangled states 
are still useful because, for example, cryptographically secure 
keys can be constructed from them fl4ll . they can enhance 
the fidelity of conclusive teleportation — hence, the telepor- 
tation power — of another bipartite state lfl5ll . they manifest 
irreversibility in the entanglement formation-distillation cy- 



cle (hence, also of interest in thermodynami cs) 11611 . and they 
have applications in channel discrimination IU7I1 . It is, there- 
fore, intriguing to see how or whether BE can appear in quan- 
tum many-body systems. Recently, it has been shown that 
a thermal environment can naturally induce the thermal tran- 
sition from FE to BE in quantum many -body systems lfl8ll . 
Additionally, in a simple XY spin chain, varying the mag- 
netic field can generate bound entangled GSs 11911 . In a recent 
experiment 12011 . an optical four-qubit state with BE has been 
realized and then characterized fully. 

There is vast literature on QPTs induced by a transition 
from separability to FE or a transition in the FE region (in 
the sense of bipartite reduced density matrices of the GSs). 
In this work, in contrast, we are interested in QPTs in which 
the transition happens because of a change in the type of en- 
tanglement, specifically, from FE to BE (or conversely), or 
from separability to BE (or conversely). This is particularly 
interesting, because one of the vastly used measures of bi- 
partite entanglement in the context of QPTs, negativity |2lll . 
is inherently unable of detecting BE, hence it would fail to 
signal a corresponding criticality. There are, however, some 
other computable (but not necessarily conclusive) insepara- 
bility tests for recognizing BE, based on the cross norm lf22tl 
and permutation 12311 criteria. In the case of quantum spin 
chains, a complete set of spin squeezing inequalities has also 
been introduced which can be used for experimentally detect- 
ing BE of many -body thermal states [24]. More appealing and 
arguably a stronger identifier than entanglement measures for 
detecting QPTs is provided b y th e ground-state fidelity (GSF) 
l25tl or its second derivative l26ll . fidelity susceptibility (FS). 
GSF and FS have been successfully used in the context of 
QPTs in correlated quantum many -body systems |27ll . 

We address the types of QPTs we are interested in through a 
reverse-construction method. Given a many -body density ma- 
trix q, we construct a spin-one chain Hamiltonian H whose 
GS has q as the reduced density matrix of specific adjacent 
sites. Thereby, by tuning the Hamiltonian parameters in some 
suitable intervals, we will have GSs representing QPTs ac- 
companied by the transition of their entanglement type from 
BE to FE or separability. The technique is in principle rela- 
tively general and is illustrated by three examples. We should 
remark that we are not aiming to construct physically real- 
izable Hamiltonians. Thus, here we are not concerned with 
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the possible appearance of many-body (three-body or higher 
order) interactions in the constructed Hamiltonians. This, al- 
though might seem to restrict the realistic applicability of the 
approach, is not a fundamental issue. There exists a power- 
ful perturbative ("gadget") machinery to construct arbitrary 
fc-body (effective) interactions from two-body Hamiltonians 
JH. 



II. CONSTRUCTION OF THE HAMILTONIAN 

Consider a closed chain of 2N identical quantum systems, 
e.g., spins, each of which has a d-dimensional Hilbert space. 
Let us assume that we are given a global state | ^) whose two- 
site reduced density matrices (TSRDMs) have some preset 
properties. For example, for odd i, Tr i i+1 [\^f)i,...,2N(^\] 
must be compatible with an already given Qi.i+i , where 
+ l = {1, . . . , 2N} - {i, i + 1} (we denote the TSRDMs 
for even i by §'). For simplicity, let us suppose that there is 
a symmetry by which all g^i+i {g^ are the same for all 
odd (even) i. That is, the matrix form of g 12 {q' 2 3) can be 
considered as a representative of all ft^+i (q[ with odd 
(even) i. We also assume that 

Qi,i+i = E^=i k\vi) l , i+ i(v l \ +0 x J2 d j=i K 

f?M+i = £*=i KWi)i,i+i{v[\ + x J2*=i K \Wj)x,i+i(Wj\, 

(1) 

where < X t < 1 for all 1, (fjiv) = (v'^v^,) — 5 n t, 

KK) = ( w j\ w 'f) = S iS> and ( v *\ w i) = ( v i\ w 'j) = °- That 
is, Eqs. ([T]) are the spectral representations of gi^+i and q\ i+1 

(including degeneracies), where the vectors {\w J )} d =l K and 

{\w' )} d J J[ K constitute their null eigenspaces. 

Our goal is to construct (or, exactly speaking, reverse con- 
struct) a corresponding Hamiltonian H for which this | "f) is a 
GS. We follow similar steps as in Ref. I29I1 . In the following, 
we restrict ourselves only to positive semidefinite Hamiltoni- 
ans. Thus, for any arbitrary |*), we have (^\H\^) > E > 
0, so that if for a |\I>), (^\H\ty) = this is an eigenvector of 
H with zero eigenvalue, namely, the GS of the Hamiltonian. 
Moreover, we are interested in nearest-neighbor two-body (or 
at most three-body) interactions. That is, the sought-after H 
has the following general form: 

H = Si odd H^i+l + Ysi even ( 2 ) 

It then can be shown that 

<*l# l*> = £« M Tr[H i!i+m , i+1 ] + £. eveD Tr[Hl !i+1 gi !i+1 }. 

(3) 

Thus, if for odd (even) i we construct their corresponding 
Hi t i + i (H^ from the null eigenvectors (or any combina- 
tion or subset of the vectors in the null eigenspaces) 
({K>}), as 

Hi, i+ i = Y^jejh l 3 ] \ w j)i,i+i( w 3\ (4) 



(and similarly for H' il+1 ), where J C {l,...,d 2 - A'} 

(J' C {1, . . . , d 2 - K'}) and hf {h'f) are arbitrary nonneg- 
ative numbers, then all the terms in the summations of Eq. (0) 
would vanish. This of course is a sufficient condition, not 
always necessary. Overall, with this choice for H, we have 
= 0, implying that |*) is aGS. 
Put briefly, then, the construction boils down to finding the 
null eigenspaces of Qi$ and ^2,3- As long as, for the given 
I 'J'), the nullities of these density matrices are nonzero, one 
can find a nearest-neighbor two-body Hamiltonian describing 
the interactions on the chain. Otherwise, the above recipe 
fails to provide coupling between all the neighboring parti- 
cles on the chain. This, however, does not necessarily im- 
ply nonexistence of other compatible two-body Hamiltonians. 
To remedy the issue with our construction, given that is 
still the state with the desired properties, a possible way is to 
consider higher order interactions, e.g., three body couplings. 

We should consider Qi.i+i,i+2, or in general ft^+i i+L (for 

some L), rather than ft^+i (for odd i). Now if we repeat 
the Hamiltonian construction recipe, we end up with L-local 
Hamiltonians iA,i+i,...,i+L obtained from the null eigenspace 
of Qn+i,...,i+L- Thus, for example, the following Hamilto- 
nian 

" 2^, //'.< I i+L; (5) 

where i could be odd or even, provides a many-body interac- 
tion for which 1^) is a GS. 

As a general example for the Hamiltonian construction, let 
us consider the following global state: 

|*)l,...,2JV = Yl!Ll VXK)l,2K)3,4 ■ ■ ■ K)2N-1,2N- (6) 

This specific form implies that except for £>j,j+i, for odd i, 
the reduced density matrices of all other pairs of sites are 
separable l30ll . The proof is straightforward, e.g., by using 
the Schmidt decomposition for the bipartite states \v t )i,i + i = 

T,e \/cFK H )^ W )2+i> for odd i. Thus, the state g x>2 fully 
captures the pairwise (i.e, two-site) entanglement property of 
the whole chain. For example, if gi 2 is bound entangled, no 
two-site FE then can be distilled from Our examples in 
Sec.|IV]are all of the type of Eq. (O; besides, they have zero- 
nullity g' 2 3 , hence H will be many bodied. 

A few remarks about the recipe for the construction of the 
Hamiltonian are in order here, (i) There could be many Hamil- 
tonians (with a set of given requirements, such as symme- 
try, range of interactions, and so on) for which a given l^) 
is a GS. Our construction provides only one class of such 
Hamiltonians. (ii) By construction, our recipe implies that 
underlying symmetries of the state g would also carry over 
into the Hamiltonian H^HU]. (iii) Note that |*) is not nec- 
essarily the unique GS of the constructed H. For exam- 
ple, each term |wj)i,2 ■ ■ • \v t )2N-i,2N (for any i) is also a GS 
(one can construct other combinations which are GS). Thus, 
when g is mixed (i.e., K > 1, which is a generic case), the 
GS is degenerate. This means that if we consider the sys- 
tem with the Hamiltonian H, at zero temperature the sys- 
tem can be everywhere in the ground eigenspace spanned by 
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{K)i,2 ■ • ■ K)2iv-i,2jv}f£i and other possible GSs. Since 
this degeneracy is symmetry-driven, we can circumvent it by 
perturbing the system with a small symmetry-breaking term, 
i.e., H' = H + eHsB, and then allowing e — > 0]. For our 
purposes here and for the sake of clarity, however, the working 
GS is assumed to be | "J) (which, for example, can be consid- 
ered to be singled out by some selection rule), (iv) In most 
physically interesting cases, we have L = 2 (2-local Hamilto- 
nians). However, since our goal in this paper is just to provide 
a framework for BE QPTs, we will not be concerned with the 
possible appearance of many-body interactions in the Hamil- 
tonians. Additionally, there exists a powerful perturbative 
gadget mechanism to construct arbitrary three-body Hamil- 
tonians from two-body interactions lr28Tl . (v) Often (but not 
necessarily always), if the state g (or depends on some 
controllable parameter (or a set of parameters), say, a, so will 
the Hamiltonian H(a). In the following, after briefly intro- 
ducing GSF as the main tool for analyzing QPTs, we shall il- 
lustrate our construction through three different examples fea- 
turing BE. 

III. ORDER PARAMETERS 

Entanglement is a rather universal means of studying QPTs 
HH. The ground-state fidelity (GSF) JH has been shown to 
be a very powerful tool that is capable of identifying various 
types of QPTs, usually independent of symmetries, system- 
dependent details, or the nature of underlying quantum cor- 
relations; hence, it can be used as a suitable "order param- 
eter." Recalling that negativity, 1 a widely used entanglement 
measure, fails to signal BE, and that other measures which 
may signal BE are usually inconclusive, GSF or its derivatives 
seem appropriate supplemental tools for our analysis. 

The rationale behind the notion of the GSF is that, at zero 
temperature, the GS describes the whole characteristics of a 
typical quantum system. Hence, when a QPT occurs at some 
point a c , a sudden change between the behavior of the GS 
slightly before (a c — A) and after (a c + A) this critical point 
may take place, which in principle can be captured by the fol- 
lowing fidelity: 

^(a c ,A) = K*(o c -A)|*(o c + A))|. (7) 

The second derivative of the GSF with respect to A, calcu- 
lated at A = 0, is called the fidelity susceptibility — hereafter 
denote d by S — and per se is a powerful measure for detecting 
QPTs [26]. For the state |*(o)) as in Eq. ©, T reads 

■F=IES'=i y/K(a-)\ t .(a + )(v t (a-)\v t ,(a+)) N \, (8) 

in which a± = a c ± A, and N (A) will be taken to be rela- 
tively large (small). 



1 Negativity of a bipartite state £) 43 is defined as N(q) = (||£> T ' 4 || l — 1 )/2, 
where q £ , = g c b t ad (in a given basis) is the partial transposition and 
= Tr[VAfA] is the trace norm HJ. 



Here we remark that the reduced fidelity (susceptibility) — 
denoting the similarity between the density matrix of a sub- 
system of the GS before and after the QPT point — can also 
be a good local order parameter for identifying symmetry- 
breaking iH and topological (H QPTs. 

There are also other, relatively less universal measures (than 
the GSF, but in some sense stronger than negativity or concur- 
rence |0]) to recognize (particularly) BE associated with the 
TSRDM of the GS, e.g., the realignment measure (shortly, 
realignment) Mr Ir22tl . It has been shown that if a bipartite 
density matrix g is separable, then, for the trace norm of its 
realignment g R (defined element-wise as g R b cd = g ac ,bd)' we 
should have || g R \\ i ^ 1. It has been shown that in some cases, 
even if M(g) = 0, the positivity of M R {g) = {\\g R \\i - l)/2 
can be a signature of BE j22Tl . Nevertheless, the nonpositiv- 
ity of Mr does not necessarily imply separability i22ll . Re- 
alignment together with negativity have been shown to set 
a lower bound for another widely used entanglement mea- 
sure, concurrence C(g) H]. More specifically, it has been 
proven that for an M\ x M2 (M\ «C M2) bipartite quantum 
state g, we have C(g) ^ 2caf, x max(M(g),MR(g)), where 

CMi = J 2 1 [Ah (Mi - 1)] II3411. Moreover, realignment has 
been recently used in devising a scheme for directly measur- 
ing entanglement of arbitrary states 1I35I1 . Recently, a crite- 
rion more powerful than realignment for iV-partite systems 
(when N is even) has been proposed ll36ll . Based on this cri- 
terion, which provides a sufficient condition for the separabil- 
ity, a sparable bipartite state g is the one for which the relation 

is satisfied. Here, g^ and gs denote the reduced density ma- 
trices associated with the subsystems of g. Since in the special 
examples we study in the next section Mr and Msr are equal 
up to the multiplicative factor 1/2, in the rest of the paper we 
make use of Mr as our main entanglement measure. 



IV. EXAMPLES 

In the following, we choose three different bipartite bound 
entangled states, and construct their corresponding Hamilto- 
nians (all spin-one chains) as in Sec. [TT] Our goal is to find 
QPTs where, for some parameter regions, a transition from 
BE to FE (or vice versa) or from BE to separability (or vice 
versa) may occur. 



A. Example I 

Consider a 3 x 3 bipartite system, consisting of two spin-1 
particles (each of which can have spins {0, ±1}, with the cor- 
responding orthonormal states {|0), |1), |I) = | — 1)}), whose 
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FIG. 1: (Color online) T of example I, where N = 10 8 and A = 
1(T 6 . 



joint density matrix g(a) has the following form l37ll : 
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where a € [0, 1]. This state, for all the parameter interval, 
has been shown to have PPT. Thus, negativity naturally fails 
to capture its entan glem ent, whereas closer inspection by the 
range criterion I37ll38ll shows BE for a <E (0,1) and separa- 
bility for a = 0, 1. 

To proceed with the construction of Sec. |IIJ we need the 
following nonzero eigenvalues and the corresponding eigen- 
vectors of g(a): 

Ai = A 2 = A 3 = A 4 = A 5 = o/(l + 8a), 
A 6 (7) = (1 + 3a T Vl - 4a + 7a 2 )/[2(l + 8a)], (10) 

and 

h) = |10), \v 2 ) = |11), |«s) = |01>, 

\v 4 ) = \01), \v 5 ) = |10>, (11) 
\v 6{7} ) = [S±(\U) + |00» T v/T=^|Il) + 7t |11)] /Vz, 

where 



6± = Vl -4a + 7a 2 ± (l - 3a), 

7± = Vl - 4a + 7a 2 ± 2a, (12) 
Z=l-a 2 + 2 < 5 2 = + 7 2 . 

Since in this example nullity of p2,3 is zero — which is also 
the case for the next examples — we increase the range of the 
interaction to three, i.e., L = 3, and calculate the null eigen- 
vectors of 0i,2.3- This density matrix has 16 null eigenvectors, 
from which only the following ones depend on a: 

\W t ) = |toi)®|i) , \U X ) - |tfl2>®|i), (13) 



Mr 
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0.0000 i 




0.0 0.2 0.4 0.6 0.8 1.0 a 

FIG. 2: (Color online) Nr(q) and W({?) (dashed line) for example I. 
Note that Mr reaches its maximum at a ~ 0.2365. 

where* G {1,0,1}, \w\) = 1 11) + .g|ll) - |11) (uptoanor- 
mahzation factor), \w 2 ) = ff '|ll) + |00)+g"|ll)-(l+3')|ll) 
(also up to normalization), g = — a)/(l + a), g' = 
-2/(3 + a), and g" = a/1 -a 2 /(3 + a). Here, K) and 
1 1^2) are the nu ll eigenvectors of gi j2 . The other vectors in the 
null eigenspace can be rewritten as {|i>i), |«2)} ® {|0), |1)}, 
{|«3>, |«4>}®{|1>, and |f6>®{|l), |0)}. Since our aim is 
to study QPTs relative to the external parameter a, the relevant 
vectors in the null eigenspace of ^1,2,3 are me ones m Eq. ( TT3T >. 
To have more symmetry, we construct the local Hamilto- 
nian only from (the normalized) \Wi) and \W\), with the 
same coupling strengths. That is, iJi.i+i^+2 = |W / i)(W r i| + 
I Wi)(Wi I (for odd i), hence H = Y,% odd ^i,i+i,i+2, or more 
explicitly as the following: 

H = E» odd [8 + 4{S* • S w ,SfSf +1 } + - 4(5, • 5 J+1 ) 2 
-JiSf + J-iSi +1 2 + ,hSf 2 S- +1 + JiS-Sf +1 2 
+9 2 (SfSf +1 2 - SfSf +1 ) + ig(Sf +1 + Sf) 
+4({5* 5»} i+ {5* + .9((Sf 2 - 5? )5? +1 2 

-^f(Sf +1 2 + ^ +1 )))]^ +2 2 , (14) 

where J 1{2) = 4(g ± 1) , J 3 ( 4) = -g(2 =p .9), and 5 = 
(S" 1 ", 5^, 5 Z ) is the spin-1 operator, and (A, B)+ = AB+BA. 

We now use the GSF T to find possible QPT points (i.e., 
criticalities) of this model (Fig. Q]). There are two points 
at which T shows nonanalyticity: a Cl = 1 and a C2 = 
1/3. Direct calculation shows that at both of these criti- 
cal points, the behavior of the GSF for even and add TV is 
different: ]im( a , A )->(i,o) H®, A ) = [8 + (-l) w ]/9 and 
lim (Q , AH(1/ 3, 0) '.F(a,A) = [9 + 2(-l) Ar ]/ll (notice the 
order of this limit taking). At a Cl the density matrix (O is 
separable; i.e., this critical point corresponds to the transition 
from BE to separability in the TSRDM of the GS of the sys- 
tem, Eq. (O, which, according to Fig. [2 can also be identified 
by AFr. Moreover, from Eq. ( fT2l . it is simple to find that at 
a Cl , 5+ = and 8- =7+; hence, \vq) and \vj) become sep- 
arable and maximally entangled states, respectively. Besides, 
here, J3 = J4 = and 3\ = — J 2 . On the other hand, g at a C2 
has BE. But, as clearly shown in Fig. [2] this criticality cannot 
be captured by AFr (or its derivatives). 

Upon closer inspection, the behavior of the GSF at very 
large (small) N (A) reveals another possibly critical point in 
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FIG. 3 : (Color online) T of example II TEq. ifToll vs a, for A = 1CT 5 
and JV = 10 9 , 2 x 10 9 , 4 x 10 9 , 6 x 10°, and 8 x 10 9 from top to 
bottom. 

this model at a C3 = 0. Since g becomes separable at this point, 
this criticality corresponds to the transition from separability 
to BE in the TSRDM of the GS. 



B. Example II 

Consider the following family of three-parameter 3x3 
states: 

X(a,(3,a) = |(1 - a - (3 - a) I <g> 1 + aP 0Q 

+ §(Pio + P 20 ) + f(Poi+Pii+f 2 i), (15) 

where P mm , = (l/ mm , ® l)|*+)($+|(C^, m ' ® 1). 

= E™»=o e2flmm 'V">(K + "i")mod3| (m,m' e 
{0,1,2}), and |$+) = ELoKO/v^. The entanglement 
properties of x(a, /3, a) have been studied, by using a geo- 
metric entanglement witness, in Ref. ll39ll . and x appears to 
have PPT BE for a ^ 0. 

To enable the method introduced earlier for the construction 
of the Hamiltonian H, it is required that the TSRDM g have 
at least one null eigenvector that depends on a (set of) param- 
eter (parameters). Note, however, that since the eigenvectors 
of x{ a i Pi a ) are parameter independent, this state cannot be 
used as g for our mechanism. Instead, we consider the follow- 
ing modification: 

g = p[qX TA + (2 !r- 1 ® 1] + (1 -P) |$ + )(* + |- (16) 

The parameters must be chosen such that the positivity of the 
state g is guaranteed. Also, by suitably choosing the param- 
eter q, we can make some of the eigenvalues of g vanish. In 
addition, introducing the parameter p is another way by which 
we can adjust the amount of entanglement of g so that a criti- 
cality occurs in our desired point. 

Because of the special form of x> me eigenvectors of g, 
apart from a dependence on the sign of pq, do not depend on 
the absolute values of p and q. This means that the parent 
Hamiltonian H has a continuous degeneracy with respect to 
the parameters p and q. In the rest of the discussion, we as- 
sume < p < 1. 




-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 a 

FIG. 4: (Color online) M and Mr of example Il-a vs a. 



There are two proper choices for the parameter q: q ± = 
2/ [2a + 2/3 - a ± 3y/{P - 2a) 2 + a 2 }, either of which can 
make three eigenvalues of g vanish. For q^, the sets of three 
null eigenvectors u ± = {v^}^ =1 are 

K ± )-(6f|10) + |01))/^/l + (6f)2, 

\vf) = (bf\ll) + |Tl))/ v /l + (6 2 ± ) 2 , (17) 

\vf) = (bt\01) + \10))/^l + (bf) 2 , 

in which bf = [aT sgn(g) a/P - 2a) 2 + a 2 ]/(2a - (3) and 
b% = bf(a —t —a). Note that when we choose q ± , the set 
{v T } corresponds to nonzero eigenvalues Xi.2,3 — P {q [a — 
2 (a + /?)]+ 2 ± 3q sgn((7)^(/3-2a) 2 +a 2 }'/18. The other 
eigenvalues of g, which we require to be nonzero (otherwise 
the degeneracy of the parent Hamiltonian would increase), are 
A 4 = {{2q{a + (3)-qa~8}p + 9}/9<md A 5 = A 6 = p[2q(a + 
0) — qa + l]/9, corresponding to the following eigenvectors: 

\vi) = (|11) + |00) + |11))/V3, 

|« 6 > = (-|11> + |ll»A/2, (18) 
| V6 ) = (-|11) + 2|00)-|11))/V6. 

Also note that, irrespective of the values of the parameters in 
Eq. (O, g' ii+1 (the TSRDM of sites i and i + 1), for even i 
has an empty null eigenspace. Thus, as explained in the pre- 
vious example, we consider pi.2,3 for the construction of the 
Hamiltonian H. The null eigenspace of pi.2,3 is spanned by 
13 vectors, from which only the following six are parameter 
dependent: 

l^i ± o) = K ± )®{|i),|o)}, 

|t^i) = l» a ± )®{|l>,|l>}, (19) 
\V ± l ) = \vi)®{\0),\l)}i 

where ± correspond to q ± . The other null eigenvectors are 
{|01), |10>> ® {|11), |Tl)} ® |0), {|10), |01)} ® |1) and 
\vs) ® |0). For odd i, we write the local Hamiltonian as 
Hi, i+W = \W±)(W±\ + IV^X^I and obtain the global 
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Hamiltonian 

H = J2i odd[ J l( S i ' + ' Si+i,SfSf +1 }+ - SfS z i+1 ) 

i T oz 2 oz 2 i T / QZ2 QZ QZ QZ 2 , QZ QZ \ 

+ J20 i D i+1 -1-^3^ D i+1 Dj -f- D i — 

+ Sf+Sf +l 2 ]{t-S z2 ) i+2 , (20) 

where 

J! = 2rfv£, J 2 = -2( M * + ^) 2 , J 3 = (/if) 2 - K*) 2 , 

Mi = fef/v/i + ^f) 2 ' = /HfS- (2D 

H represents a three-parameter family of Hamiltonians 
whose parameters are determined so that the positivity of g 
is guaranteed. In the following, we study the entanglement 
properties of the GS of this Hamiltonian in the vicinity of its 
critical points, for two special cases. 



1. Example H-a 

When a = (6 - 6)/21, (3 = -26/21, and a = (5 - 
2b)/7 [i.e., a = (1 + a)/6 and /3 = (-5 + 7a) /m 
X(a, /3, a) is reduced to another bound entangled state 14011 
X = [2|$+)($+| + 6(7+ + (5 - 6)<7_]/7, in which <7+ = 
(| 10) <10| + |01)(01| + |ll)(ll|)/3, and cr_ is the swap of a + . 
For 6 € (3,4], \ h as BE. Let us assume that q = q + and 
p = 3/4. The parent Hamiltonian then will be in Eq. < [2Qb 
with superscript +. The GS (O is written with the set ir and 
its corresponding eigenvalues, {Ai, A2, A3 } | g=9 + , plus eigen- 
values {A4, A5, Ag}| 9=g + and eigenvectors of Eqs. < fT~8b . In 
this way, the GSF (O can be obtained readily. Figure [3] shows 
the finite-size behavior of T vs a, for A = 10~ 5 and different 
N. As is seen, T has a minimum at a c = for all values of N. 
When increases, T decreases; however, the rate of decreas- 
ing at a c = (ViV) is maximum. This behavior can be a sign 
of criticality in the system. At this point, bf = bf and nega- 
tivities of \vf ), \vf ), and \vf) attain their maximum values. 
On the other hand, as seen in Fig. |H a c = is in the region 
where Af(g) vanishes; hence, Af fails to signal criticality in 
the system. From the positivity of A/r(o) around a c , it is evi- 
dent that a = is the point where the bipartite entanglement 
of the GS experiences transition from BE to separability — Mr 



detects this critical point. Furthermore, despite the nonanalyt- 
icity of the first derivative of Mr at points g goes from the FE 
to the BE region, there is no critical behavior in T at these 
points. This implies that the transition from BE to FE cannot 
be responsible for QPT in this system. 



2. Example H-b 

Suppose that a = (—1 + 3a)/6, ft = (1 + 3a)/7, q = q~, 
andp = 0.76027256. According to our previous explanations, 
the Hamiltonian H (the GS | \I/ ) ) is constructed from the set v~ 
[v + and its corresponding eigenvalues {Ai, A2, A3 } |q =(? - , to- 
gether with the eigenvalues {A4, A5, A6}| q=9 - , and the eigen- 
vectors (1181) 1. Thus, the global Hamiltonian of this model will 
be H in Eq. ( f20b with superscript — . In Fig. [51 S/N has been 
plotted in terms of a for different values of N . There exist two 
points at which this system may exhibit criticality. At a Cl = 
5/6, S/N becomes nonanalytic. Indeed, for arbitrary val- 
ues of p at this point, we have liiri( a ^)-s.(5/6,o) F(a, A,p) = 
{11 + 10[-1 + (-l) N ]p}/n. At a Cl , since 2a = f3, then 
^1(2) — * °°- Therefore, vectors {v^} in Eq. ( TTTb at a Cl are 
completely factorized. Furthermore, one should note that at 
this point, in the global Hamiltonian H, /iw 2 ) ^ 1 an< ^ tnen 
u t(2) S oes ^° zero - ^ s clearly shown in Fig. |6l this proba- 
ble QPT occurs when g is completely free entangled, and of 
course negativity does not show it (note that this result is due 
to the existence of a trace norm in the definition of negativ- 
ity since the first derivative of one of the eigenvalues of g TA 
has a discontinuity at rj Cl and therefore can signal it). Another 
possible QPT point of this system is a C2 w 0.229650. As the 
behavior of S/N in Fig. [5jshows, the rate of the drop in T at 
this point is maximum, which may signal a QPT in the system. 
The interesting fact is that this critical point corresponds to a 
change in the type of entanglement of the TSRDM of the GS. 
In fact, since a C2 is the initial point of the range at which M is 
zero and Mr is positive (see Fig. |6j, it represents a transition 
from FE to BE in the two-site states. As Fig. [6] shows, the 
first derivative of Mr at a C2 is discontinuous, implying that it 
detects this criticality. 
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FIG. 7: (Color online) M and Mr of example III vs a. 



FIG. 8: (Color online) T of example III vs a, for N = 10001 and 
A = 0.01. 



C. Example III 

As the final example, consider the following 3x3 density 
matrix: 
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in which 

g = l + a + a 2 , b = ^a + g + 1 , 7 = (1 + g) 2 , 
w = 2.a6, cr = fe 2 (.o-l), /i= (l + a)(l + 3.g + a),(23) 
v = 2gb 2 , P = 2^6(1 + a) , r\ = -o 2 a, e = (2g + a) 2 . 

After some straightforward algebra, we find the following 
nonzero eigenvalues and their corresponding eigenvectors: 

Ai = 3/10, A 2 = ••• = A 8 = 1/10, 
K (2) ) = [|ll)±fe|10) + (l + o)|ll)]/(l + a), 
\v 3 ) = [ a ( a + l)|ll)+6 2 |00)-a|ll)]/6 2 , (24) 
\ Vi ) = |11), |« B ) - |01), K) = |0T), 
\v r ) = |11>, |«b) = |10>. 

For the same reason as in the previous examples, i.e., the nul- 
lity of 02,3 constructed from Eq. (O is zero, we work with 
Qi. 2. 3 to construct the Hamiltonian. The null eigenspace of 
Pi. 2, 3 is spanned by 

\W t ) = [-(1 + o)|ll) + o|00) + |11)] g) | t ), 
W) = [-(^)|H) + [00) + i|ll)] ® |0), (25) 
|t0 = (a'|ll) + |ll))®|0> > 

together with {|l0),|ll)} ® {|1),|0)}, {|01>, |01)} ® 
{|1),|I)}, |11) O {|0),|I)}, and |100), where 1 G {1,1} 
and a 1 = 1/(1 + a). We consider the local Hamiltonian as 



ilj ! j_l_i i j_l_2 = \U)(U\, for odd i, hence 

H = J2i oddl^il 2 + i s i ■ Si+i,s^Sf +1 }+ - (St ■ s l+1 ) 2 
-sf - st + 2 + {s*, sy} i+ {s*, sv} i+1+ ] 

+J 2 (Sf 2 S- +1 2 + S*S* +1 ) + J 3 (Sf 2 Sf +1 + S*Sf +1 2 )} 

x(l-S* 2 ) i+2 , (26) 

in which 

Ji = -4a' a", J 2 = (a' + a") 2 , J 3 = (2a' 2 - 1), (27) 

and a" = y/\ - a' 2 . 

Figure |7] shows the behavior of negativity M and realign- 
ment Mr vs a. At a £ (—2, —0.4], M remains zero, whereas 
Mr is negative (with a nonanalyticity at a c = —1). As ex- 
plained before, although for a separable state Mr is nonpos- 
itive, the converse is not necessarily true. We could not con- 
clude whether any BE can exist in this specific region of a 
because we do not know of any other measure stronger than 
Msr (whose behavior — as explained in Sec. [Ill] — is the same 
as Mr up to the multiplicative factor 1/2) for detecting BE. 
We, nonetheless, show that a criticality occurs in this region 
exactly at the nonanalyticity of Mr. Figure [8] depicts the GSF 
(O of this model, showing a QPT at a c = — 1. At this point, 
li m (a,A)-s-(-i,o) J- (a, A) = [2(— 1) N + 3]/5, and nonanalyt- 
icity at a c appears for odd A^s. 

Finally, we remark that we investigated the behavior of the 
single- and two-site fidelity (susceptibility) of our examples 
around their critical points. The results show that in most of 
the cases, the total behavior of these local order parameters 
are similar to that of their global counterparts, hence they are 
also capable of detecting quantum critical points associated 
with our models. For example, the single-site fidelity in ex- 
ample I shows a sharp drop at a Cl and the two-site fidelity 
(susceptibility) in example Il-a (example Il-b) can identify all 
the QPT points. Nevertheless, there exist some critical points, 
e.g., a C2 in example I, for which these local measures do not 
herald quantum criticality. 

V. SUMMARY AND DISCUSSION 

Here we have investigated the role of bound entanglement 
in quantum phase transitions. To this aim, we have presented 
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a method to construct quantum spin-chain Hamiltonians in 
which a quantum phase transition can be accompanied by a 
change in the type of entanglement (of the two-site reduced 
density matrices of the corresponding ground states) from 
bound to free or separability. This method per se is fairly 
general and can be applied to various scenarios. Given a form 
for the desired two-site reduced density matrices (with their 
engineered entanglement properties), by utilizing their spec- 
tral properties, we have outlined a reverse construction for a 
compatible pure state for the whole chain. Additionally, we 
have suggested a method to obtain a class of Hamiltonians 
for which the constructed pure state is a ground state. At 
zero temperature, then, by varying the Hamiltonian param- 
eters (carried over from the very parameter dependence of 
the given reduced density matrices), a desired quantum phase 
transition at some engineered points can take place. To be able 
to detect possibly relevant criticalities (due to the underlying 
entanglement-type changes), we have used the ground-state 
fidelity and the realignment criterion. 

One direction for extending our results is to delve into a 



possible connection between our construction of Hamiltoni- 
ans and the methods of quantum inverse scattering |4T1i . Al- 
though these methods have clear distinctions, there might be 
a rich common framework through which one can theoreti- 
cally engineer quantum models with prescribed entanglement 
properties for ground states. 

Overall, our study may hopefully spur further interest in 
bound entanglement and the role it may play in quantum phase 
transitions or other areas in which a potential application or 
advantage could be anticipated. For example, quantum com- 
munications 1 1711421 14311 and controllable generation of entan- 
glement in engineered physical systems 14411 could be possible 
beneficiaries. 
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